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Abstract 



A theory of a confined two dimensional electrolyte is presented. The 
positive and negative ions, interacting by a 1/r potential, are constrained to 
move on an interface separating two solvents with dielectric constants ei and 
€2- It is shown that the Debye-Hiickel type of theory predicts that the this 
2d Coulomb fluid should undergo a phase separation into a coexisting liquid 
(high density) and gas (low density) phases. We argue, however, that the 
formation of polymer-like chains of alternating positive and negative ions can 
prevent this phase transition from taking place. 

PACS numbers:64.70.-p, 05.70.Fh, 64.60.-i 

I. INTRODUCTION 

Over the last decade the need to understand coulombic criticality has provided a new 
impetus to the study of electrolyte solutions. The current wave of exploration can be traced 
to the pioneering experiments of K. S. Pitzer et. al, who have first reported a surprising 
finding that Coulomb interactions might belong the mean-field universality class This 
suggestion has not gone unchallenged and, in fact, later experiments are consistent with a 
crossover from mean- field to Ising universality class very close to the critical point |J. The 
crossover, if indeed it exists, is much closer to the critical point than for any other known 
fluid. 

From the theoretical point of view it is very hard to justify anything but Ising critical- 
ity [[IU . The goal for the theorists must then lie in a seemingly simpler task of finding why 
the crossover region for electrolytes is so narrow. Unfortunately even this programme has 
failed to produce any satisfactory explanation |J . Most calculations suggest that the critical 
region for electrolytes should be comparable to that of a Lennard- Jones fluid. To further 
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confound the mystery, Monte Carlo simulations are once again pointing in the direction of 
mean-field criticality ||. 

A seemingly unrelated problem concerns the disappearance of the anticipated liquid- 
gas transition in a system of dipolar hard spheres (DHS). Since the DHS is the simplest 
realization of a polar fluid, for a long time it has been believed that it must exhibit a 
liquid-gas phase separation. It came, therefore, as quite a surprise when the Monte Carlo 
simulations failed to located this singularity |7[. Instead what they found was that as the 
temperature was lowered, the dipolar particles aligned, forming polymer-like chains. Since 
these chains interact weakly ||, it has been argued that the liquid-gas transition must be 
driven entirely by the free unassociated dipoles 0. In fact a critical density of free dipoles 
is needed for the phase transition to occur. At low temperatures, where most of the theories 
localize the transition, formation of dipolar chains depletes the density of free dipoles bellow 
the critical threshold necessary for the transition to take place ||. 

In this paper we shall study a third model which we hope might span the gap between the 
three dimensional Coulomb gas and the dipolar hard spheres, and thus shed some additional 
light on the criticality in these interesting and important systems. Our new model consists 



of a neutral electrolyte confined to a two dimensional plane flQfl . This can be visualized 
as oppositely charged surfactant molecules adsorbed to a water-oil interface. An example 
of such a system is cetyltrimethyl amonium-hydroxy naphthalene carboxylate (CTAHNC), 
which is composed of two surface active parts, CTA + and HNC~ fllT] . 

We shall argue that unlike the 3d electrolyte [|12| , the confined 2d plasma might not phase 
separate. Instead as the temperature is lowered, chains composed of alternating positive and 

negative ions will begin to form (... H 1 — ...) Q. Just as dipolar chains, these new clusters 

interact weakly between themselves. However, they diminish the concentration of free ions 
bellow the critical value necessary for the transition to take place. 



II. THE MODEL 

Our system consists of an ionic fluid of N + positive and N_ negative ions confined to a 
plane of area A located at z — 0, separating two different solvents occupying the half-spaces 
at z > and z < 0. We shall restrict our attention to the neutral electrolytes, for which 
N + = iV_ = N/2. The solvents are treated as uniform mediums with dielectric constants 
62 and e% for z > and z < 0, respectively. The ions of both species are modeled as hard 
spheres of diameter a and charge ±q located at the center. It will be convenient to define 
the dimensionless particle density as p* ± = p±a 2 , where pi = Ni/A. The total density of ions 
is p = p + + p-. 

All the relevant thermodynamic information is contained in the free energy density, 
/ = F j A. Unfortunately due to the complexity of interactions, no exact expression for / 
can be found. We shall, therefore, attempt to construct the approximate free energy using 
the most relevant contributions. These can be divided into entropic and electrostatic, 

f(T, p+, p_) = fent(T } P+,P-) + fel(T, P+, /£>_) . (1) 

The entropic (mixing) free energy is given by 

/3/i2(T, P+ , p_) = [p + ln( P+ A 2 ) -p + + p- ln(p_A 2 ) - p_] (2) 
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where A = J 2n% 2 /inksT is the thermal de Broglie wavelength of the ions. The second term 
in Eq. ([D is due to electrostatic interactions. It is important to note that the electrostatic free 
energy is purely correlational, since the mean- field contribution is zero. To calculate f e i, let 
us fix one ion at the origin. Adopting the cylindrical coordinate system (g, (p, z), the central 
ion is located at g = 0, z = 0. Due to the electrostatic interactions the other particles will 
arrange themselves within the plane in accordance with the Boltzmann distribution. Since 
no charge is present in the regions z > and z < 0, the electrostatic potential there must 
satisfy the Laplace equation. Appealing to the azimuthal symmetry and taking into account 



the fact that the potential should vanish at infinity we find [13 



</> 1 (g,z) = A 1 (k)J (kg)e kz dk for z<0 (3) 
Jo 

and 

roo 

(f) 2 (g,z)= A 2 (k)J (kg)e- kz dk for z>0 , (4) 
Jo 

where Jq(x) is the Bessel function of order zero. 

The functions Ai(k) and A 2 (k) are determined by the boundary conditions: continuity 
of electrostatic potential, 

Mq,o) = Mq,o) , (5) 

and discontinuity of displacement field across the z = plane, 

(€ 2 B 2 (g,z) - e^^g^)) ■ n = Arra eff (g) , (6) 

where a e ff(g) is the surface charge density and n is a unit vector normal to the interface, 
pointing from region 1 to 2. The continuity of electrostatic potential results in Ai(k) = 
A 2 {k) = A(k), while Eq. @ requires that 



2 / i, A(k)J (kg)dk = 47Ta ^ /(g) , (7) 



o 

where D — {e\ + e 2 )/2. The surface charge density |T3| is given by 

veff(e) = a s (g) + . (8) 
lug 

The term o s is the charge density of the "ionic cloud" around the central ion, 

o s (Q)=qp + e~ (3q<t, -qp-e + ^. (9) 

Eqs. (0) and (H) are only valid in the limit a — * 0. Unfortunately it is highly nontrivial to 
take into account the boundary condition associate with the finite ionic size. To circumvent 
this difficulty, we shall first solve all the equations in the point particle limit. Then to 
account for the finite particle size, we shall replace the bare charge of the central ion q, in 
Eq. (P), by an effective charge Q, q — > Q. The effective charge Q will be determined by the 
condition of an overall charge neutrality, 
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27T / a s (g)gdg = -q. (10) 

J a 

For the 3d electrolyte this procedure leads to an electrostatic potential identical with the 
one found from the exact solution of the Debye-Hiickel equations with the appropriate hard 



core boundary conditions [14]. In the present geometry the procedure outlined above will 
only be an approximation, although we believe a rather good one. 

In the spirit of Debye-Hiickel theory |L4| we shall now linearize the Boltzmann factor in 



Eq. (g). The surface charge density then reduces to 

271 

where k s = 2np* /T*a is the inverse Gouy-Chapman length and T* = ksTDa/q 2 is the 
reduced temperature. Eq. (0) can now be solved yielding the expression for electrostatic 
potential, 

4>{q,z) = % r -J—j o (k )e- k ^dk. (12) 

L) Jo k + k s 



For z = 0, the in plane potential agrees with the one obtained by Velazquez and Blum flO 
and can be conveniently rewritten as 

<PAq) = — — , (13) 

The subscript " >" is included to stress that for rigid particles this form will be appropriate 
only outside the hard core exclusion region, g > a. The function t v (x) is defined as |TD 



t v {x) = 1 7i-[H v (x) - N u (x)} (14) 



where H v (x) and N v (x), are the Struve and the Bessel functions of order u, respectively [p~5 



The charge neutrality, Eq. ([T0|), together with Eqs. ( PD and (p~3|) , determines the effective 
charge Q, 

Q = (15) 

K s aTi{K s a) 

In the limit of a — ► 0, the renormalized charge reduces to the bare charge q. 

To obtain the electrostatic free energy, we require the potential inside the excluded vol- 
ume region, 4><{o)- Unfortunately the procedure presented above leaves this undetermined. 
The only statement that we can make is that 0<(f?) must be of the form 

Mq) = ^- q +C(q). (16) 

For a three dimensional unconstrained electrolyte, C is a constant. This, however, is not 
the case in the present geometry and C(g) is a function of position. In particular the value 
C(0) is the potential that the central ion feels due to the presence of other particles. In their 
earlier study, Velasquez and Blum [Kj approximated C(g) by a constant which they then 
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determined by requiring continuity of electrostatic potential across the exclusion boundary 
4><{a) = </) > (a). This, however, is a very rough approximation since there is nothing to 
prevent C(g) from being a very strongly varying function of position. To avoid this difficulty 
we shall use an alternate method of obtaining the value of C(0). To this end we note that 
the potential at the center of a circularly symmetric charge distribution is 

977- roc 

C(0;p) = — I a s ( Q )d Q , (17) 

where in order to emphasize that the potential depends on the density of ions, we have 
explicitly included p in its definition. Using Eq. ([□]) with the potential given by Eq. (|13D 
we find 

C(0;p)= q r^dz . (18) 
Dan(k s a) JK S a z 

The excess chemical potential can now be calculated straight forwardly by appealing to the 
Guntelberg charging process fL6| . We find /i± = qC(0; p)/2. The chemical potential of the 
positive and negative ions is 

= £ =Mf) + «- (19) 

The critical point is determined by the conditions, dp±/dp = and d 2 p±/dp 2 = 0, which 
reduce to 

2T*r 2 + In - Ir - Tiro = 
ItI - AI t qT i - 3t t? + Irl + 2r 2 n + x 2 Ir\ + In + x 2 rf + r 2 = (20) 

with 

/(a;) = f°° I^ldz , (21) 

where Solving Eqs. ( p0|) , the critical point is located at T* = 0.0517386 and 

p* c+ = p *_ = 0.00121988. 

III. DEBYE-HUCKEL-BJERRUM APPROXIMATION. 

Clearly the low temperature at which the critical point is located should make us worry 
about the approximations which have been adopted. Certainly at such extreme conditions 
the linearization of the Boltzmann factor, Eq. (0), is no longer valid. Fortunately, all is not 
lost. Evidently, linearization of the Boltzmann factor in Eq. ([J) diminishes the weight of 
configurations in which the oppositely charged ions are in a close proximity. It is possible, 
therefore, to correct for the omitted non-linearities by explicitly allowing for the formation 
of "clusters". These clusters are assumed to be in a chemical equilibrium with the free 



unassociated ions, and their density is determined by the law of mass action [17,18]. The 
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most basic such cluster is a dipole formed by a (H — ) pair. Within the simplest version of this 
theory the dipoles are treated as non-interacting ideal specie. For 3d electrolyte this Debye- 
Huckel-Bjerrum approximation (DHBj) ||,|J has proven extremely successful, predicting the 
location of the critical point in close agreement with the Monte Carlo simulations [|12[ . 

The total density of ions p can then be subdivided into that of free unassociated 
monopoles p\ = p + + p_, and of dipolar pairs P2, with p = pi + 2p 2 . In the spirit of 
DHBj theory we shall first treat the dipoles as ideal non-interacting entities whose concen- 
tration is determined by the law of mass-action, p + + fi_ = p 2 . The Helmholtz free energy 
density is then given by 



PL 



DHBj 



fentiPl, P2,T) +f el (p 1 ,T) 



where the entropic contribution is 



Pfent(Pl,P2,T) 



P2 In 



P2 



Pi In 



>iA 2 



Pi 



(22) 



(23) 



The £2 is the internal partition function of a dipolar pair, 



Z 2 (T;R) = 2ir [ e 

J a 



gdg 



(24) 



where fiUi = a/T*g is the electrostatic potential between the associated ions. In order to 
evaluate £ 2 (T;i?) we must define the distance R at which two ions can be considered to 



be associated. Following Bjerrum []T7 we choose the value of Rsj at which £ 2 (T;i?), as a 
function of R, has an inflection point, Rsj = a/T* With this, the integral in Eq. (|24"D 
can be evaluated to yield 
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na 2 e b 



b s e~ b (Ei(b) - Ei(l) + 2e) - 6(1 + b) 



(25) 



where b = 1/T* and Ei(x) is the exponential integral function. 

Since within the DHBj approximation the dipoles are treated as ideal, the electrostatic 
free energy only depends on the density of free monopoles, p\. Substituting the free energy 
into the law of mass-action, the density of dipoles is given by 



P*2 = &P+P-e 



(26) 



The charge neutrality requires that p + = p_ = pi/2. The expression for excess chemical 
potential was calculated in Section 2, p:± = j3qC(0; pi)/2, and the inverse Gouy-Chapman 
length is now k s = 2np\/T*a. 

Within the DHBj approximation the dipoles are treated as ideal specie, therefore, they 
cannot influence the location of the critical point. Thus, the critical temperature must still 
be T* = 0.0517386 while the critical density of monopoles must remain p\ c = 0.00243975. 
Substituting these values into Eq. (|26"D, we find that the density of dipoles is p* 2c = 1.08847, 
which is extremely high. If there are so many dipoles is it also not possible that there will 
be higher order clusters as well? 
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IV. LINEAR IONIC CHAINS 



Unfortunately as soon as we get to clusters of three ions the calculations get extremely 
difficult. The basic problem is the internal partition function of the higher order clusters, 
which can no longer be evaluated exactly. Furthermore, while it is evident that for a cluster 

of three ions the low temperature configurations are chain-like, (H h) or ( — I — ), this is 

far from obvious for a neutral cluster of two positive and two negative ions. The entropy 

favors a chain-like configuration (H 1 — ), while the energy favors a compact square cluster. 

Which will gain in the critical region is hard to say. We note, however, that exactly the same 
situation was encountered for dipolar hard spheres PJT9|]. In that case, in the vicinity of the 



critical point, the chain like configurations dominated. Since it is much easier to study the 
linear clusters they, therefore, will provide a starting point for our analysis. 

We begin by supposing that at low temperatures our system will be composed of 
monopoles of density p\ and chains of n monomers with densities p n . Once again in the 
spirit of Bjerrum we shall first treat the clusters as ideal non-interacting species. The particle 
conservation requires that 

oo 

p=Y / np n . (27) 

n=l 

Consider an alternating chain composed of t positive and s = n — t negative ions. The 
partition function for such a cluster is, 

^-waLU^'^ "- (28) 

where fl n is the configurational volume and U n is the total energy of interaction between 
the ions forming a chain, 

n 

U n (r 1: r 2 , r n ) = ^ Vijfai) • ( 29 ) 

i<j 

At zero temperature chains are rigid and the particles are in contact with one another. The 
displacement vector between two ions i and j is ■ = a\i — j\x, where x is the unit vector 
along the chain. The electrostatic energy of interaction between the ions of the chain can 
be evaluated exactly yielding 

K = $-S{n) , (30) 



Da 



where 



s(B) = £(z!£jp*). (31 ) 

fc=i K 

For nonzero but small temperatures, deviations from • occur. These fluctuations can be 
taken into account by making a Taylor expansion of U n around the ground state up to 
quadrupolar order 
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J2 n 



£4 = t/„°-^£(-i) i+J { 



('I') 2 



+ 



9 If 13 



|f?.|3 



}• 



(32) 



Here the subscripts x and ?/ represent the components along the chain's direction and per- 
pendicular to it, respectively. Since we are assuming small deviations from the ground state, 
for each ion i we shall consider only fluctuation contributions to £ n arising from the inter- 
actions between the nearest and the second nearest neighbors. Choosing as the unit bases 
x and y, vectors parallel and perpendicular to the direction of the chain in the plane z — 0, 
the nearest and the second nearest displacement vectors can be written as 



(0 



ail + A» 



COS( 



sin i 



(33) 



and 



Rf = a 



1 + Aj) cos <pi + (1 + Ai+i) cos 

(1 + Aj) sin0j + (1 + Ai+i) sin0 i+ i 



(34) 



respectively. Here A is the radial and is the angular deviation from the relative positions 
in the ground state. 

Substituting Eqs. (^) and ( p4|) into Eq. (0), for small fluctuations the electrostatic 
energy becomes 



U„ 



1 



n-l 



n-2 



knT T 



k=\ 



Afc + Afc + i — ^(0/c 



(35) 



In the low temperature limit, configurational integral, Eq. (p8|), can be performed explicitly 
yielding, 



,2n-2 ' 
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2 2n 7r f r ^-2 e -5(n)/T* /or n > 3 . 



(36) 



The prefactor is the result of thermal fluctuations while the exponential is due to the ground 
state energy. 

The condition of chemical equilibrium between the monopoles and the n-chains is ex- 
pressed through the law of mass action \i n = The chemical potential for monopoles 
is given by Eq. (|T9| ) with p — > p 1; while for n-chains the chemical potential is, 



[i n = k B T In 
The law of mass-action reduces to 



2n 



Pn 



(37) 



(38) 



At the level of approximation that we have adopted, the chains are treated as non- 
interacting ideal species. This means that just like in the case of dipoles in Section 3, they 
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cannot affect location of the critical point. Therefore T* = 0.0517386 and p\ c = 0.00243975. 
Substituting these values into Eq. (|38|) we find that the sum in Eq. ( pT|) diverges. In fact 
according to the Cauchy criterion, the sum in Eq. ([27]) will converge absolutely if and only 
if, 



A = lim p* n 1/n < 1 . 



(39) 



Using Eq. (BSD, 



A = 2p*v / vrT* 3 exp 



a + I(x)/t 1 (x) 
2T* 



(40) 



where a = 2 In 2 and x = n s (pi)a. Inserting the critical parameters into the expression 
above, we find that at criticality A c = 5.82706, and the Cauchy criterion is strongly violated. 
Therefore, the critical density of monopoles lies outside the radius of convergence of Eq. (p7|). 
This means that for any finite density p, the density of monopoles never reaches the threshold 
necessary for the phase separation to occur, pi < p lc . 



V. RESULTS AND DISCUSSION 

We have presented an argument which suggests that a confined 2c? electrolyte should not 
phase separate. Instead, just as for the case of dipolar hard spheres, as the temperature 
is lowered, the ions will associate forming polymer-like chains of alternating positive and 
negative monomers. Clearly our argument is based on a number of assumptions. First, we 
have supposed that in the critical region the linear chains predominate over the compact 
clusters. This is far from obvious. If the compact clusters have lower free energy than the 
chains, they can provide nuclei for the condensation and the gas-liquid phase separation. 
A second assumption implicit in our calculations is that the chains interact only weakly. 
This is somewhat easier to justify. Consider an infinitely long rigid line of alternating 

charges (... H 1 1 — ...). Suppose that the center to center distance between the nearest 

neighbors is a. It is then possible to show that the potential produced by such a line of 
charge decays exponentially, ip(r) ~ exp(— irr/a), where r is the distance perpendicular 
to the chain. Thus, the interactions between long polymer-like clusters should, indeed, 
be quite weak. However, the shorter chains can still interact sufficiently strongly to drive 
phase transition [M. Finally, even if the formation of chains prevents the liquid-gas phase 
separation, it does not forestall other kinds of phase transitions from taking place ||21|| . At 
sufficiently high densities, the two component plasma will crystallize exhibiting a pseudo long 
range order. At low temperatures and densities, where the polymer-like chains predominate, 



the Y-like defect formation |22[ can lead to a coexistence between two phases, one with high 
and the other with low concentration of defects |[23|| . 

After this work was completed and submitted for publication, the referee drew our at- 
tention to a recently published simulation by Weis, Levesque, and Caillol (WLC) |24j] of a 
2d ionic fluid. Indeed, these authors found coexistence between high and low density phases. 
WLC estimated the critical temperature to be T c « 0.04, which should be compared with 
our Debye-Hiickel prediction of T c = 0.052. Furthermore, in the low density phase, WLC 
found predominance of chain and ring-like clusters. The high density phase resembles a 
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percolating gel-like cluster J24|. Although there is a phase coexistence, it is difficult to asso- 
ciate it with a traditional liquid-vapor transition. Instead the coexistence resembles more a 
sol-gel transition in polymer systems. The task for theorists must now be to quantitatively 
understand the phase transition found by WLC. We hope that the current paper will provide 
a first step in this direction. 
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